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We consider the quantization of a Friedmann-Robertson- Walker universe. We derive a reduced 
square-root Hamiltonian by choosing the scale factor as time variable and quantize the theory using 
Pauli matrices d la Dirac. From the resulting spinor equation we show that there is no semiclassical 
wave packet that avoids the big bang singularity. Our work raises the question concerning the 
relationship between the choice of time and singularity avoidance. 
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I. INTRODUCTION 

To understand space and time at the Planck scale we 
have to quantize general relativity (GR). An attempt to 
achieve this goal consists in the program of canonical 
quantum gravity. Here, one has to capture the dynamics 
of a globally hyperbolic spacetime in terms of a Hamil- 
tonian formulation of GR. Following the work of Pirani 
and others [l], 0], Bergmann and others Jj], as well as 
Dirac [4j, Arnowitt, Deser, and Misner [J wrote down 
the sought-after formulation in terms of metric variables 
in 1962. Using this result, the quantization can be per- 
formed according to Dirac's work on constrained Hamil- 
tonian systems [6J, which gives rise to what is called 
the Wheeler-DeWitt equation 0, [1] . Another method is 
referred to as reduced phase space quantization. This 
procedure has been extensively developed in work by 
Kuchaf [I, [Io| in the 1970s. The idea is to solve the 
Hamiltonian and diffeomorphism constraints on the clas- 
sical level. This brings them into a form that, after a 
Dirac quantization, gives rise to a functional Schrodinger 
equation, called a Tomonaga-Schwinger equation (Tll[l3|. 
Finally, the program of loop quantum gravity (LQG) pro- 
vides another way to quantize GR. Here, instead of us- 
ing metric variables, one works with so called Ashtekar 
variables (l3| . The origin of this approach lies in work 
performed by Sen [14| in 1982. 

There are a number of problems that come along with 
both Dirac and reduced phase space quantization. The 
most prominent is the one of time [HI, [l6|, 0]. When 
compared to the Schrodinger equation, the Wheeler- 
DeWitt equation does not have the time derivative on 
the right-hand side. Among other things, this implies 
the lack of an obvious Hilbert space to work with. Re- 
garding the reduced phase space quantization, the prob- 
lem of time seems to be more severe. First of all, the 
method fails when applied to full GR [HI]. Secondly, 
when it comes to using the technique for simple models, 
different classical time choices lead to different reduced 
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Hamiltonians. Therefore the Tomonega-Schwinger equa- 
tions associated with them are not the same, which raises 
the question about the validity of one choice of time over 
any other. Also, as the reduced Hamiltonians are square- 
root functions and in general time dependent, deriving 
and solving the functional Schrodinger equation is very 
difficult or even impossible. 

In spite of everything, the interest in the ideas and con- 
cepts mentioned above is unwaning. This is because the 
focus has been shifted from the full theory to its applica- 
tions in simple models. Most prominent are Bianchi mod- 
els, whereof the isotropic Friedmann-Robertson- Walker 
(FRW) cosmologies play the most important role. The 
reason for this is twofold. On the one hand, spatially 
flat FRW cosmologies provide us with the best model 
we have for our Universe; on the other, they state that 
space and time started with a big bang singularity. This 
is precisely where Einstein's theory of gravity is supposed 
to break down and where quantum effects are expected 
to become crucial. In fact, the expectation and hope is 
that they lead to an avoidance of the big bang. This 
has been shown to be true not only for FRW cosmolo- 
gies [ill HH H3, [23[ but also in other cases such as 
spherically symmetric spacetimes 0, HH HE H3 • How- 
ever, the matter used is as simple as ours, namely a mass- 
less scalar field. There is not a single case known to the 
author where more realistic matter has been successfully 
incorporated. 

In this paper we investigate the problem of the singu- 
larity avoidance in a flat FRW universe with a nonnega- 
tive cosmological constant and a massless scalar field as 
the source. In Sec. |TT]we state the well known classical 
theory of the model and, in particular, the Hamiltonian 
formulation thereof. It is in this section where we estab- 
lish a basis for the quantization. To be more precise, we 
solve the Hamiltonian constraint by choosing the scale 
factor of the FRW metric as time variable. This gives 
rise to a reduced Hamiltonian, which is a square-root 
function of time and the momentum canonically conju- 
gated to the massless scalar field. In Sec. IIIII we first 
perform the quantization with the help of a trick that 
Dirac used to derive his famous equation from the dis- 
persion relation in special relativity. Then we solve the 
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resulting Schrodinger equation; we construct wave pack- 
ets that are semiclassical at late time and show that they 
do not avoid the big bang singularity. Finally in Sec. 
IIVI we give a general discussion and make a comparison 
to results from loop quantum cosmology (LQC), where a 
singularity avoidance is achieved. 



II. CLASSICAL THEORY 

The model we are investigating can be characterized 
as follows: First, we use a spatially flat FRW metric 

ds 2 = -N(t) 2 dt 2 + a(i) 2 <%dxW 

with a nonnegative lapse function N and a scale factor a. 
Second, the source of our spacetime is a massless scalar 
field 0, minimally coupled to the gravitational field. We 
choose it to be a function of t only so that the symmetries 
of the metric are respected. Finally, we have a nonnega- 
tive cosmological constant A. 



A. Lagrangian formulation 

The field equations and their solutions can be derived 
from a Lagrangian formulation. One of its main ingredi- 
ents is the Ricci scalar 



R = 6 



N 2 c 



( a \ Na 
\Na) ~ N^a 



If, in addition to the lapse function, we restrict the scale 
factor to be nonnegative as well, drop an overall factor 
V representing what is called the fiducial volume and set 
k = 16irG, we can write 
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(1) 



for the Lagrangian. On the matter side, this gives rise to 
an energy density and a pressure given by 
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respectively. As far as the field equations are concerned 
we can write down two independent ones for the three 
unknowns N, a, and <p so that a choice has to be made. 
Since the velocity N is absent in the Lagrangian, it is 
most natural to determine the lapse function. We fix a 
time variable, a clock. If we set it equal to 1 then we are 
dealing with the proper time r and, according to [28| . 
the solutions to the field equations are 



a(r) = (V3Ar) , 

<Xt) = </) + In (VXr) 
V3ft ^ ' 



(2) 



for A = 0, and 
a(r) 



sinh v 3 At 



1/3 



(r) - 0o + 



/3ft 



: arctanh 



(3) 



for A > 0. The constant A = Ka 6 <fi 2 /4 is positive and <t> Q is 
specified by the initial data. The solutions we presented 
just now imply a big bang singularity for all A > (this 
holds true even for A < 0). It does not vanish for other 
choices of the lapse function as we can simply rescale the 
time axis. In other words, GR predicts a singularity for 
the model we study, independent of the lapse function 
and the cosmological constant. 



B. Hamiltonian formulation 

We are interested in studying a quantized version of the 
given model. The path we choose to get there leads over 
a Hamiltonian formulation. We initiate it with the speci- 
fication of the Arnowitt-Deser-Misner (ADM) action [f|. 
If we parametrize the metric hij on spatial hypersurfaces 
according to hij = a 2 Sij, where a is any sufficiently dif- 
ferentiate function of i, and if we Legendre transform 
the matter part of the Lagrangian in Eq. (fTJ), the ADM 
action in canonical form is given by 
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The momenta are 



Pa 



(p a a + - H) dt. 
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the nonvanishing Poisson brackets {p a , a} = 1, 
1, and the Hamiltonian 



H = N 
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is constrained to vanish. We include the fiducial volume 
V in order to keep track of all units. On the Lagrangian 
level, the vanishing of H is equivalent to the Friedmann 
equation. However, on the Hamiltonian level it implies a 
static universe and to get dynamics, a time variable has 
to be chosen. We focus on a clock defined by y/Xt — a(t) 
which, suppressing the argument p^, gives rise to the 
reduced Hamiltonian 

1/2 

I • (5) 
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Since we want the choice for the time variable to be pre- 
served as t changes, we have to have 1 = {H, A^ 1 / 2 a} 
which, using Eq. ([5]), implies 
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Had we decided to go with the negative square-root of 
the reduced Hamiltonian, we would refer to a contract- 
ing phase of the Universe (with respect to any coordinate 
time that is monotonically increasing). This can be seen 
from the definition of p a in Eq. ((4]). We stick with the 
positive square-root since the clock we defined describes 
an expanding phase. Now, for t > 0, Hamilton's equa- 
tions of motion are solved by 
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for A = 0, and 



<t>{t) = <t>o 
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for A > 0. In both cases, p<j>(t) = P4, is a positive con- 
stant. If it is given by 2V{\/ k) 1 / 2 we can recover Eqs. @ 
and ([3]) by plugging in the proper time variable r, which 
can be constructed from the lapse function. As can be 
seen from Eq. ([3]), this means that if A > 0, the expan- 
sion of the Universe is accelerating at late time. Finally, 
the equations 



R(t) = 4A 
-T° (t) 



1 



2V 2 X 3 t 6 ' 
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and T\{t) 



1 



2U 2 A 3 i 6 



imply that we have the same big bang singularity as in 
the Lagrangian formulation. 



III. QUANTUM THEORY 

Because of the square-root, the Hamiltonian in Eq. ([5]) 
is troublesome. Even though we can insert a Schrodinger 
representation for the momentum p^ and, at least for 
A > 0, directly get a self-adjoint and positive square-root 
Hamiltonian 



h a (t) = 



48U 2 A 3 A . 12 K 
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on L 2 (M.), we cannot solve the time dependent Schroding- 
er equation. As has been carefully pointed out in [29(, 
an attempt to circumvent this problem by solving the 
Klein-Gordon equation fails because of the explicit time 
dependence of the Hamiltonian. All we can do with the 
square-root Hamiltonian is to solve the time independent 
Schrodinger equation 

h a m E (t) it, <j>) = E{t)^ m (t, <f>). 

The reason for this is that, here, we can square on both 
sides. It leads to the eigenstate 

i> E{ t)(t,4>) = i>_{t)e- iV W h + 4> + (ty Vm/ \ (9) 



where the momentum V(t) is given by 



V{t) = ( T2 E{tft 2 



4U 2 A 3 A 
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and the functions i^±(t) are determined by the initial 
data. Note that the expression for V(t) introduces a 
lower bound on E(t), specified by the cosmological con- 
stant. Now, because of 



h a (ti), h a (t) 



0. 



(11) 



the eigenstate in Eq. is an eigenstate at a later time 
ti > t as well. Therefore, if the Universe is in an eigen- 
state at ti, its energy E(t{) is given by 



E(ti) 2 tf = E(t) 2 t 2 



48U 2 A 3 A 



itf-t 6 ). (12) 



This is a relation between the vacuum energies E(t{) and 
E(t). The difference between them comes from the en- 
ergy introduced by the cosmological constant. Using this 
relation and Eq. (fTI]) , we can evolve the eigenstate in Eq. 
© to t t . For t > and A > this yields 



where the integral in the exponent is 
E{t{)ti E(t)t 



(13) 



E{t\) dt[ = ., 
t 



'3k 



V(t) In 



y/U/tfPjt) + E(tj)tl 
Uj y/l2/KV(t)+E(t)t 



The case A = is discussed in [29| . This is all we can say 
about the solutions to the time independent Schrodinger 
equation for the square-root Hamiltonian. Again, be- 
cause of the square-root, the time dependent Schrodinger 
equation cannot be solved. Fortunately this is not the 
end of the story as we can quantize the reduced Hamil- 
tonian in another way. 



A. Dirac's trick 

Motivated by Dirac's analysis of the Klein-Gordon 
equation, which led him to derive the Dirac equation, 
we start with an Ansatz for the Hamiltonian operator of 
the form 



where A and B have to be Hermitian matrices. Note 
that we set t > 0. If we square the Ansatz on both sides 
and make a comparison with the squared Hamiltonian 
from Eq. ©, we get a system of equations for A and 
B that is solved by any pair of distinct Pauli matrices 



1/2 
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and transformations by an arbitrary nonsingular matrix 
thereof. Here, we set 



.4 



1 

1 



and B 



1 

-1 



The motivation for this choice will become clear further 
down. Using a Schrodinger representation for the 
Hamiltonian is now 



h a {t) = 



at 2 



at 2 



where a and (3 are chosen in accordance with Eq. (|14p . 
It is now clear that A < is not possible as it pre- 
vents the Hamiltonian from being Hermitian. If we define 
H = L 2 (R) © L 2 (R) as the Hilbert space, with the inner 
product of two states fi £ Tt given by 



we make h a (t) a sclf-adjoint operator on H. Also, the 
induced norm of a state iff = (ip, lu) 1 £ Ti is conserved. 
This can be seen from the continuity equation 



dt 



M 2 )+d [/?/i(|vf-M 2 )] =o, 



whose origin lies in the time dependent Schrodinger equa- 
tion 



h a {t)f{t,(f>) = ihdt*(t,<f>). 



(15) 



Note that the components ip and u> decouple in this equa- 
tion if we neglect the cosmological constant. The same 
decoupling occurs when t approaches zero. This is the 
limit in which the big bang singularity can be found on 
the classical level and it motivates our choice for the ma- 
trices A and B. 

The question remains whether the singularity can be 
avoided when we use the quantum mechanical formula- 
tion we built just now. To answer it we first have to 
solve the time dependent Schrodinger equation. Then we 
have to construct semiclassical wave packets and check 
whether the expectation values of the Ricciscalar, the 
energy density, and the pressure, with respect to these 
states, remain finite when t goes to zero. In the appendix 
we consider other choices for the matrices A and B and 
we show that the final conclusions remain the same. 

Before we go into the discussion on the singular- 
ity avoidance, we take a look at the time independent 
Schrodinger equation 

K{t)^ E{t) {t,<t>)=E{t)^ E(t) {t,<f>). 
The solutions are 



*B(t)(«. 
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(1 \ ■ E(t)t < 

JJ^ + (tK^ (i6) 



for A = 0, and 



Lin 3 r>M I ^+(*) e 



iV(t)<p/h 



for A > 0. As before, the functions ip±(t) are specified 
by the initial data and the momentum V(t) is defined in 
Eq. (p~0|) . Here, V{t) introduces both an upper bound on 
the negative energies and a lower bound on the positive 
energies by means of the cosmological constant. Using 
the same reasoning as for Eqs. (fl~2]l and (fl~3|) . if A = 
the energy at f ; > t is given by E{t{) = E(t)t/ti and the 
eigenstate in Eq. (fT6|) by 



**!(*,) (W) 



iE{t)t/h 



However, if A > 0, we cannot write down similar expres- 
sions for the energy and the eigenstate at ti because we 
now have 



h a (ti),h a (t) 



This implies that h a (ti) and h a (t) do not have common 
eigenstates, a property which is used for the computation 
of E(ti). Also, it means that the time ordering symbol in 
the Dyson series, which defines the evolution of the eigen- 
states, cannot be neglected. This being said, we conclude 
the discussion of the time independent Schrodinger equa- 
tion. 



B. Singularity avoidance? 

a. Case A = 0. We need to solve the time depen- 
dent Schrodinger equation. This is a task that can be 
simplified if we replace t and <fi by the dimensionless co- 
ordinates y/Xt and <p/j3 respectively. Namely, it allows us 
to rewrite Eq. (|T5)) and get 



td t + 8^ 
tdt - d 4 



*M) = o, 



(17) 



for which the solution is given by 
4>(t,<l>) =/(0-ln(t)) and u(t, 



g(4> + Ht)). 



The functions / and g are in i 2 (R), but otherwise arbi- 
trary. Therefore, we can write down a wave packet 



(0-gn-ln(t)Vf 



-0o-ltt(t)) 



with normalization constant m 
be seen from 



tt-^W- 1 / 2 . As can 

h 2 

(P0>* = hk and (p 2 )* = h 2 k 2 + , (18) 



5 



among the constants k and W the former describes the 
momentum of the wave packet, whereas the latter char- 
acterizes the quantum fluctuations of the squared mo- 
mentum operator. Since the expectation values 



and 



(</>)* = + hi(i) and ((f> 2 )^ = (<p + ln(t)) 2 + 



2 



imply that \1/ is always peaked over the classical solution 
from Eq. ([6]), with standard deviation W/V2, we have a 
semiclassical wave packet if we require k 3> 10 ~ 34 (the 
numerical value of K) and 1 > > 10 -34 (in units 
where (3 = 1). Now, the point is that this holds true 
for all t > 0. The evolution has no impact on the form 
of the wave packet and tells us only where in the (t, <j>)- 
plane iff is most likely to be found. In Fig. [T] we plot the 
probability density function |^| 2 for small and decreasing 
values of t. It indicates that the evolution is undisturbed, 
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FIG. 1: We chose </><> = 1, ko = 4, and W = 0.5. The peak 
traces the classical solution and evolves freely. 

and that the wave packet follows the classical solution. 

As far as the big bang singularity is concerned, the 
equation 



lim |*fi,< 

t-»o 



= 



implies that the probability density function vanishes if 
<fi G M is fixed and t approaches zero. However, this 
only reflects the fact that the wave packet bypasses finite 
values of (j> and follows the classical, diverging solution 
given in Eq. ©. Indeed, from Eq. (fT5|) we get 



2V 2 X 3 t 6 

i (i>lh 



(19) 
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2V 2 X 3 t 6 
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(20) 



The big bang singularity is not avoided. 

The above conclusion holds true for arbitrary wave 
packets. Any 'J/ that solves Eq. (fT7|) can only have an 
upper component if it is to be semiclassical. Since this 
component is a function of <fi — <f> — ln(i), an integration 
over <f> can be replaced by one over <p — <p — ln(i), thereby 
removing the ^-dependence. Therefore, the expectation 
value (p^)* is independent of t. Since the squared mo- 
mentum operator has a spectrum on L 2 (K) that is strictly 
positive (except for the irrelevant constant functions), the 
big bang singularity cannot be avoided. 

b. Case A > 0. Here, we work with the dimen- 
sionless coordinates (a/(6/i)) 1 / 3 t and 4>/f3 so that the 
Schrodinger equation becomes 



td t + 
i6t 3 



i6t 3 
td t - dd 



= 0. 



Its solution is given by 

i>(t, 0) = it 3 ' 2 J [j(fc)J„(2t 3 ) + y(k)Y l/ (2t 3 )]e lk ^ dk 
and 

w(t t <f>) = t 3 ' 2 f [j{k)J v+1 (2t 3 ) + y(k)Y v+1 {2t 3 )]e ik * dk, 



where v = — 1/2 — ik/3. The functions J„ and Y v are the 
Bessel functions of the first and the second kind of order 
v respectively. See [30| and [3l[ for useful details related 
to these functions. If we set 

■-^(fe-fco) 2 -ifc(0 o + i ln(A/A)) 



j(k) = ncosh ^ — kj e 



and 



y(k) = msinh y^k j e 



-(fc-feo) 2 p -ifc(0 o + i ln(A/A)) 



where n — 2 1 / 2 7r 1 / 4 IF 4 / 2 is a normalization constant, 
we get the wave packet 



<j>) = nt 



= »+ 3 / 2 



cosh 



iJ„{2t 3 ) 
Ju+i(2t 3 ) 

-Y„(2t 3 ) 
iY„ +1 (2t 3 ) 

^(fe-fe o ) 2 el fc(</»-0o-|ln(A/A)) dfc 
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sinh ( —k 
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(21) 



The choice for the functions j and y is motivated by the 
fact that when t is large, 'J/ is semiclassical. Namely, in 
this limit we have 



i cos(4i 3 ) 
"" sin(4i 3 ) 



(*-0o-£ ln(A/A)) 
2W' Z 



iko(<t>— <t>o — | ln(A/A)J {22) 
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Because of our choice of n, this wave packet is normalized 
and we can compute the expectation values 



(<f>h 



Hi 



(23) 



and 
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A comparison with Eq. shows that, at late time, \& is 
peaked over the classical solution. Also, Eq. l[2"2"|) gives 



(p</>)* ~ and 



2 1.2 



2IF 2 



(24) 



so that as in the previous case, if we require k a 10~ 34 
(the numerical value of ft) and 1 3> W 3> 10~ (in 
units where /3 = 1), we have a semiclassical wave packet. 
However, this does not hold true for all £ > 0, the reason 
being that \& has two components and if £ <C 1, they 
move in different directions with respect to the 0-axis. 
This becomes evident if we expand Eq. (f2"Tj) about £ = 0, 
from which we get 



V(t,0) 



and 



w(t,0) 



cosh r(-z/) 

. e -^i(fc-fc o ) 2 el fc(0~0 o -ln(t)) & 



-i^ J sinh^fc) T(y + 1) 

. e --^(fc-fc o ) 2 el fc(*-0o+ln(t)) dfc 



where T denotes the gamma function. A complete dis- 
cussion of this function can be found in [H E3 • Now, 
|^| 2 is a Gaussian function in <f> — <f> — ln(£), and \u>\ 2 is 
one in <f) — o + ln(£). Therefore, when £ gets closer to zero 
their peaks approach (f> = — oo and (f> — oo, respectively. 
As can be seen from Eq. ([7|), the component ij) follows 
the classical solution. However, the component u) does 
not and the wave packet is no longer semiclassical. This 
is illustrated in Fig. [21 where we plot the full wave packet 
from Eq. (|21[) for several values of £. 

The result in Eq. (|2"3")l is only true for large £. The 
expectation values in Eq. (f24j) however, can be computed 
using Eq. (|2ip and are exact for all £. Therefore, the Ricci 
scalar, the energy density and the pressure diverge as £ 
approaches zero, as it happens in Eqs. ([iU]) and (|20p. In 
conclusion, the big bang singularity is not avoided. 

The value of the cosmological constant has no effect on 
the conclusion concerning the singularity avoidance. It 
simply serves the purpose of giving rise to a time depen- 
dent mass that vanishes when £ goes to zero. The choice 
of the wave packet does not matter either. Namely, all 

have two components, ip and w, whereof for small £ 
the former has to be a function of cj) — 4> — ln(£) and the 




FIG. 2: We chose A = A, <f>o = 1, ko = 4, and W = 0.5. 
Starting at <j> = —6, we computed |$| 2 after every A(j> = 0.25. 
The wave packet is no longer semiclassical for small t. 



latter one of 4>— 4> Q + ln(£). Thus, they give rise to an ex- 
pectation value (p^)* that is time independent if £ <C 1. 
As it is also strictly positive on L 2 (R) © L 2 (R), the big 
bang cannot be avoided. 



IV. DISCUSSION 

In the present paper we quantized a flat FRW universe 
with a nonnegative cosmological constant and a massless 
scalar field as the source. To do so we first chose the 
scale factor as time variable, which gave rise to a reduced, 
time dependent square-root Hamiltonian. Then, in order 
to be able to perform the quantization, we referred to a 
famous trick that Dirac used to derive the equation that 
is named after him. This led to a Schrodinger equation 
that is a system of two coupled linear first-order par- 
tial differential equations. We solved it and constructed 
wave packets that are semiclassical at late time. Next 
we evolved the wave packets toward the point, where the 
big bang singularity can be found on the classical level. 
In case of a vanishing cosmological constant, this led to 
a quantum big bang as the wave packet remains semi- 
classical for all time. As far as positive cosmological con- 
stants are concerned, the big bang singularity could not 
be avoided either. This happened, even though the wave 
packet was no longer semiclassical in the vicinity of the 
singular point. Finally, we showed that all semiclassical 
wave packets lead to a quantum big bang. 

The clock we chose is expressed in terms of metric vari- 
ables. It is a good clock as it is closely related to the 
Hubble time, which is the most natural choice for FRW 
cosmologies. Blyth and Isham [29| considered the same 
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time variable we do. However, their spacetime has a neg- 
ative curvature and the cosmological constant is set equal 
to zero. Even though they find a reduced time dependent 
square-root Hamiltonian, and write down the associated 
spectrum and eigenstates, no statement is made regard- 
ing singularity avoidance. 

Clocks are made of matter so that it makes sense to 
try to use the massless scalar field as a time variable. 
This is what Blyth and Isham [2^| do. However, for this 
calculation they modify our model by letting the space- 
time have a nonnegative curvature and neglecting the 
cosmological constant. They compute the energy spec- 
trum and the eigenstates from the squared, time indepen- 
dent Schrodinger equation and show that no wave packet 
avoids the big bang singularity. 

Another result is found by Ashtekar and his collabo- 
rators [H [H HH, HI- Their model is a FRW universe 
with a nonnegative curvature, an arbitrary cosmological 
constant and a massless scalar field as the source. The 
quantization is performed using techniques from LQC. To 
be more precise, they use a polymer representation [32| 
for the gravitational part and a Schrodinger representa- 
tion for the matter part of the Hamiltonian constraint. 
The scalar field is identified as a time variable and a finite 
difference equation is derived, which leads to a resolution 
of the big bang singularity. 

It is really a combination of both the clock and the 
polymer representation that is responsible for the result 
above. On the one hand, the clock gives rise to a time in- 
dependent Hamiltonian; on the other, the polymer repre- 
sentation introduces an upper bound on the gravitational 
kinetic term. This causes the matter part to be bounded 
too, which leads to the singularity avoidance. 

The choice of the clock is very important in the context 
of singularity avoidance. In the LQC case [i~9ll20l. l2ll |22| . 
the scalar field is chosen as time variable. It gives rise 
to a time independent Hamiltonian and an avoidance of 
the big bang singularity is found. In our work the scale 
factor is chosen as clock, the Hamiltonian is time depen- 
dent and the singularity is frozen into the model. (We 
think that a polymer representation of the matter part of 
the Hamiltonian does not resolve the singularity because 
replacing p^ by sva{fj.p^) / n in Eq. (fT4|) does not remove 
the divergence in the small t limit. However, this poly- 
merization is a conjecture which requires further study 



to establish.) This also occurs in [29] and in [33|, where 
the Gowdy universe is considered. In each case, the sin- 
gularity is unavoidable as soon as the clock is introduced. 

What clocks can we choose if the scalar field is massive? 
Except for the simplest cases, any clock leads to a time 
dependent Hamiltonian. Also, at least in the context 
of reduced phase space quantization, the Hamiltonians 
are square-root functions, and quantizing them is very 
difficult. Finally, as far as solving the Wheeler-DeWitt 
equation is concerned, it is not obvious which clocks can 
be chosen through the WKB approximation [34| . 

Altogether, the ultimate question is and remains what 
can be concluded from minisuperspace models. 



APPENDIX: MATRIX REPRESENTATIONS 

First of all, we have to note that the Schrodinger equa- 
tion (fT5T> can be written as 

(-ifirfdu, + a'e 3t ) *(t, 4>) = 0. 

Here we factored out the matrix A and replaced t and 
<j> by the dimensionlcss coordinates ln(-\/Ai) and <j>/(3 re- 
spectively, set a' = a I A 3 / 2 , do = dt and d\ = d$, then 
used the gamma matrices 



What we wrote down is the Dirac equation in two space- 
time dimensions. It follows from its covariance that the 
results we derived are independent of the choice we made 
for the matrices A and B. 
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